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For quantum systems described by finite matrices, linear and affine maps of matri- 

> ' 

ces are shown to provide equivalent descriptions of evolution of density matrices for 



(N 



a subsystem caused by unitary Hamiltonian evolution in a larger system; an affine 
map can be replaced by a linear map, and a linear map can be replaced by an affine 



' map. There may be significant advantage in using an affine map. The linear map 

cn ' 

cO ■ is generally not completely positive, but the linear part of an equivalent affine map 

(N 

can be chosen to be completely positive and related in the simplest possible way to 
the unitary Hamiltonian evolution in the larger system. 



We are accustomed to the use of linear maps of matrices to describe evolution of 
density matrices in the dynamics of open quantum systems 

was pointed out recently 10| that affine maps might be used as well. I will show here that 
the linear and affine methods of description are equivalent for systems described by finite 
matrices: an affine map can be replaced by a linear map, and a linear map can be replaced 
by an affine map. There may be significant advantage in using an affine map. A linear map 
that describes evolution of density matrices for a subsystem caused by unitary Hamiltonian 
evolution in a larger system is generally not completely positive [11], but the linear part of 
an equivalent affine map can be chosen to be completely positive and related in the simplest 
possible way to the unitary Hamiltonian evolution in the larger system. The equivalence 
demonstrated here is for a finite-dimensional Hilbert space. It may be that one kind of map 
will work where the other does not when the Hilbert space is infinite-dimensional. There 
too, it may be advantageous to have two alternatives. 
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Consider a quantum system described by NxN matrices. Let L be a linear map of 
NxN matrices to NxN matrices; it takes each NxN matrix Q to an NxN matrix L(Q). 
Let K be an JVxiV matrix. The map M that takes each NxN matrix Q to the matrix 

M{Q) = L{Q)+K (1) 

is called affine. Let p and a be NxN density matrices and let 

r = qp+(l-q)a (2) 

with < q < 1. Then r is a density matrix and 

M(r) = gL(p) + (1 - g)L(a) + K 

= q[L(p) + K} + (l- q )[L(a) + K] 

= qM(p) + (l-q)M(a). (3) 

If M(p) and M(a) are density matrices, then M(r) is a density matrix and it is related to 
M(p) and M(a) the same as it would be if M were linear. The property of linear maps that 
has physical meaning for density matrices is provided by affine maps as well. 

In the iV 2 -dimensional linear space of NxN matrices we can find a basis of N 2 Hermitian 
matrices F M for p — 0, 1, ... N 2 -l such that F is 1 and 

Tr[F M F„] = N5^. (4) 

If p is a density matrix, then 

N 2 -l 

1 / 

P 



(5) 



with (F a ) = Tr[F a p], and 



1 / N 1 \ 

M(p) = ^(F a )L(F a )J+K 

1 / ^ \ 

= -l L (l) + NK+ ^(F Q )L(F a )J 



(6) 



Compare this with the result of a linear map that takes each NxN matrix Q to an NxN 
matrix Q'. It gives 

N 2 -l 



'=17(1'+ Ew.) 
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The affine map and the linear map give the same result for all density matrices if 

l' = L(l) + NK, F' a = L(F a ). (8) 

Specification of a linear map means independent specification of its action on each basis 
matrix 1 and F a for a = 1, 2, .... N 2 -l. We can choose 1' and F' a to match any affine map. 
Conversely, we can choose L(l), K, and L(F a ) to match any linear map, but the choice of 
L(l) and K is not unique. 

Specifically, suppose the density matrices are for a subsystem of a larger system. 



Evolution of these density matrices caused by unitary Hamiltonian evo 



ution in the larger 
This linear map 



system can almost always be described by a linear map of matrices jll|. 
is generally not completely positive, but the linear part L of the equivalent affine map is 
completely positive when L(l) and K are chosen so that L(l) is 1. We can see this from the 
paper of Jordan, Shaji and Sudarshan if L(l) is 1, then L is their linear map with zeros 
for the parameters that involve mean values of quantities for the larger system; and when 
these parameters are all zero, their linear map is completely positive. For the two-qubit 
example that they work out in detail, if £(1) is 1 then L is the linear map with a\ and a,2 
both zero; they observe explicitly that if a± and a 2 are zero, the map is completely positive 
for all t. In general, if L(l) is 1 then L is the linear map with the parameters d^ all zero. 
For any initial state of the subsystem, there is a density matrix for an initial state of the 
larger system that gives zeros for the parameters d^ and is a product of the density matrix 
for the subsystem and a density matrix for the other part of the larger system; all that is 
required is zeros for mean values of quantities for the other part of the larger system. When 
the parameters d^ are zero, the map can be obtained from unitary Hamiltonian evolution 
starting from an initial product state for the larger system. This implies that the map is 
completely positive. In fact, if L(l) is 1, then for each matrix p for the subsystem, density 
matrix or not, 

L{p) = ^Tr R {e~ im p®l R e im ] (9) 

where e~ tHt .. e tHt gives the unitary Hamiltonian evolution of density matrices in the larger 
system, R denotes the other part of the larger system (the remainder or rest of the larger 
system, which could be a reservoir), so Ir is the unit matrix for R and Tr R is the trace over 
the states of R, and M is the dimension of the Hilbert space for R; this holds for all the basis 
matrices except 1 regardless of how £(1) and K are chosen, and it holds for 1 when L(l) is 
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1. Thus L is related in the simplest possible way to the unitary Hamiltonian evolution in 
the larger system. 
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